Impurity-induced frustration in correlated oxides 
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Using the example of Zn-doped La2 Cu04 , we demonstrate that a spinless impurity doped into a non-frustrated 
antiferromagnet can induce substantial frustrating interactions among the spins surrounding it. This counterin- 
tuitive result is the key to resolving discrepancies between experimental data and earlier theories. Analytic 
and quantum Monte Carlo studies of the impurity-induced frustration are in a close accord with each other and 
experiments. The mechanism proposed here should be common to other correlated oxides as well. 
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Impurities are known to be an effective tool to locally per- 
turb quantum systems, thereby revealing important informa- 
tion about their microscopic interactions and correlations flJJ] . 
A well studied example of a strongly correlated quantum sys- 
tem in which effects of such impurity doping can be investi- 
gated is La2CuC>4 — one of the most important cuprate super- 
conductor parent compounds. In its pristine form, this ma- 
terial is a two-dimensional (2D) spin-i Heisenberg antiferro- 
magnet (AF) [2]. It is believed that the substitution of Cu 2+ 
(S = |) ions by spinless Zn 2+ represents a good realiza- 
tion of the site-diluted Heisenberg hamiltonian 0, El Si • 
In this Letter, we demonstrate that there exists a signifi- 
cant qualitative correction to the dilution picture. Impuri- 
ties can induce substantial frustrating interactions between 
nearby spins. Not only does this effect explain discrepancies 
between experimental data and the dilution-only theories for 
La2Cui_2;Zn K 04, but it may also be important for a variety of 
other phenomena in diluted magnets and doped Mott insula- 
tors. Our mechanism for such an effect should be common to 
many charge-transfer insulators, including oxides of transition 
metals. 

We propose that the presence of extra degrees of freedom 
due to oxygen orbitals necessarily results in frustrating terms 
in the corresponding low-energy spin hamiltonian of the Zn- 
doped system, which are absent in the dilution-only models. 
Utilizing quantum Monte Carlo (QMC) and analytic T-matrix 
approaches, we calculate the doping dependence of the stag- 
gered magnetization for such a low-energy model. We show 
that this model, with the parameters appropriate for the CuC<2 
planes given by a three-band Hubbard model calculation, nat- 
urally explains experimental data. 

Experiments and theories. — Comprehensive studies of the 
problem of La2Cu04 diluted by spinless Zn impurities have 
been performed using neutron scattering, magnetometry, and 
NMR (NQR) on the experimental side [3, 4], and QMC and 
T-matrix approaches of the diluted Heisenberg model on the 
theoretical side S |a 01- These studies allow for extensive 
cross-checks. The unbiased QMC data agree with the T- 
matrix results closely up to x ~ 15%, s_upporting the validity 
of the latter in the low-doping regime ||5l |6(] . However, there 
are serious discrepancies between theoretical and experimen- 



tal results. Fig.Q]shows the average magnetic moment M per 
Cu site versus the Zn doping fraction x. The experimental 
data are always below the theoretical curves. The slope 
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at small x represents the rate at which the order parameter 
M is suppressed by individual impurities due to enhanced 
quantum fluctuations. The inset of Fig. [TJ shows a large 
discrepancy — a factor of approximately two — between the 
theoretical and experimental results. This indicates that the 
dilution-only theory significantly underestimates the impact 
of the impurity on the quantum spin background. 

One might attempt to explain the disagreement by suggest- 
ing that longer-range ( J2, J3, etc.), and ring-exchange interac- 
tions should be included in the model for the undoped CUO2 
plane. Such terms are generally present in the low-energy spin 
models derived from the Hubbard model [9], and they do lead 
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FIG. 1: (Color online). The on-site magnetic moment M(x) per Cu 
normalized by its undoped value vs Zn doping x. Symbols show 
neutron scattering |3] (squares) and NQR |4] (diamonds and circles) 
data, with dashed lines their best fits. The dotted line is the classical 
CS— > 00) result. Solid lines are the QMC [5] and the T-matrix results 
@], respectively. Inset: slope function (see text) vs x @] 
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FIG. 2: (Color online) AF bonds in the (a) dilution-only model, (b) 
model with impurity-induced frustration (two out of six frustrating 
bonds are shown), (c) The CuC>2 plane with a Zn impurity, (d) 
Schematic local level diagram in the impurity-doped system. 

to a reduction of M, However, since the order is suppressed 
already in the undoped system, this mechanism is unlikely to 
enhance fluctuations specifically due to dilution. Using an ex- 
pansion of M in the dilution fraction x and in the extended 
interaction J 2, one obtains at small x and J2: 



R(J 2 ) « R{0) [l + A 



J2 
J 



(2) 



where A < 1 ifioll and i?(0) is the theory slope from Fig. Q] 
Thus, a large correction to the slope of the M(x) /M(0) curve 
in the extended model is only possible if J2 ~ J, which is be- 
yond the realistic range in the cuprates where J2/J is at most 
of the order of 10% lllll . A recent study |Tl'] has shown that 
while extended interactions are important for explaining the 
lower absolute value of the staggered magnetization, they are 
not able to explain the large initial slope in the M(x)/M(0) 
dependence. Thus, one must seek another explanation. 

Extra interactions. — The dilution picture seems natural for 
modeling the replacement of a magnetic Cu site by a mag- 
netically inert Zn; see Fig. Ufa). However, for the dilution- 
only picture to be valid, the Zn-site must remain electroni- 
cally inert at energy scales up to the Hubbard U. Since it is 
the three-band Hubbard model that describes the real CuC<2 
plane and other transition-metal oxides [12], the states on the 
oxygen orbitals also become important. They hybridize to Zn 
and remain involved in virtual hoppings between surrounding 
Cu-sites, see Fig. |2fc), facilitating extra superexchange cou- 
plings that connect further neighbor Cu-sites. Thus, the spin- 
less impurity, in effect, leads to a cage of frustrating interac- 
tions around itself, with four J' Zn and two J Zn ; see Fig.[2jb). 
Qualitatively, the impurity-doped system is not equivalent to 
the site-diluted Hubbard model with electronically inert im- 
purity sites, but rather to the t-e-U model, 

a visual representation of which is given in Fig. |2fd). The 



t-U part is the usual Hubbard model, which at half-filling re- 
duces to the Heisenberg model at t 2 /U order. The higher- 
order terms are negligible (~ t 4 /U 3 ) if t <C U, The dis- 
tinct physics is brought into play by the model OJ when 
the energy cost at the impurity site e is less than the Hub- 
bard gap. In that case, virtual transitions through the impu- 
rity level will cause superexchange interactions of order of 
~i 4 /e 3 . Taking e = U/2 and U/t = 10 for an estimate leads 
to J' Zn /J ~ (t/U) 2 {U/e) 3 ~ 0.1. The total impact of the 
impurity-induced frustrating interactions per impurity is then 
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• 0.6 J, which, as we will show below, is 



enough to explain the discrepancy between dilution-only the- 
ory and experiments. As is discussed above, the correspond- 
ing Hubbard terms of the 4th order are smaller and do not 
disturb the order parameter specifically due to dilution. 

For the realistic values of the CuC<2 plane parameters, map- 
ping of the three-band Hubbard model to the single-band one 
can be done using the cell-perturbation approach [13] which 
does not require the smallness of the Cu-O hopping t p d with 
respect to the charge-transfer gap A lfl4ll . In this approach, 
locally hybridized states on Cu and surrounding O's are diag- 
onalized exactly and the three-band model becomes a "multi- 
orbital" Hubbard model with the effective "Cu" states con- 
nected by effective hoppings. Since the lowest states in the 
multi-orbital Hubbard model are the same as in the single- 
band one (i.e. the lowest two-hole state is the Zhang-Rice- 
like singlet) the equivalence of the two models can be justi- 
fied lfl3l1 . In this approach, even if Zn is inert electronically, 
the remaining O-like states in the Zn-C<4 cluster can facilitate 
couplings between neighboring Cu spins, Fig.|2c). 

We extend this approach to the Zn-doped case and per- 
form a detailed microscopic calculations of J' Zn and J Zn 
Il5tl . First, we fix the parameters of the three-band model so 
they yield the experimental value of the Cu-Cu superexchange 
J~0.13eV [13]. Since the electronic parameters of Zn states 
are not known precisely |T^, 17], we vary them substantially 
as shown in Fig. [3] for a representative set of the three-band 
model parameters. Our Figs. Oa),(b) show how the energy 
Ezn of the lowest effective "Zn" state depends on the energy 
of the bare Zn-level, Az n , and the hybridization, tzn-o, re- 
spectively. The effective Hubbard energy U e f / is also shown 
to demonstrate the validity of the qualitative level structure in 
Fig. |2d) and to support our model (O. Figs. [3a),(b) show 
that the electronic levels of Zn and hybridization with them 
are important in lowering ez n and enhancing J' Zn and J Zn . 

Our Figs.[3c),(d) show the A^„ and tzn-O dependence of 
the total impurity-induced frustrating interactions per impu- 
rity J' tot / J. Given the uncertainty in Az n (from 2-3eV, 111 611 . 



to 5eV, 11711 '). the total frustrating effect can be estimated to 
be between J' tot ~ (0.2 — 1.0) J. Individual J' Zn and J Zn are 
in the range of 3-15% of J. Counterintuitively, the interac- 
tion across the impurity ( J Zn ) is greater than the next-nearest 
neighbor interaction (J' Zn ) due to partial cancellation of the 
super-exchange and the ring-like exchange involving Zn and 
three Cu sites on a nearest-neighbor plaquette, see Fig. |2fb). 
This results in a stronger bond between copper spins across 
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tion in lattices with N = L x L sites: 
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FIG. 3: (Color online) (a) and (b): energy ez n of the lowest state on 
the effective "Zn"-site (Zn04 cluster) vs Azn and tzn-o, respec- 
tively, (c) and (d): J' tot = 4J^ n +2J^ n vs Azn (for several values of 
tzn-o) and tz-n-o (for several values of Az n ), respectively. Inset 
(d): J'zn/J'zn vs t Z n~o- In (a-d), t pd = 1.5eV, A = 3eV, [7 d >A. 



the Zn-site, with the ratio J Zn / J'z n — 2 ^ 4 in a wide range 
of the three-band model parameters, see inset in Fig. [3d). Al- 
together, the three-band model provides support to our idea 
and gives an order-of-magnitude estimate of the parameters. 

Low-energy model. — With this microscopic insight, one 
should model the Zn-doped CuC>2 plane by the Heisenberg 
model with random impurities not only causing more fluctua- 
tions by cutting links [3], but also connecting the nearby spins 
in a frustrated way. Thus, the effective model is: 



(<i> (v)' 



J 



5> 

m" 



S,- (4) 



where the first term is the dilution-only model with 5=1/2 
spins for all the sites except the impurity sites (where 5 = 0) 
and the sums over the J' Zn and J Zn bonds are taken around 
the impurity sites only, as shown in Fig. [2jb). 

Suppression of the order parameter. — We have investigated 
the model © by means of the analytical T-matrix as well as 
unbiased QMC techniques. The former is based on the di- 
agrammatic treatment of the corresponding linear spin-wave 
theory with exact calculation of the scattering amplitudes off 
the impurities and subsequent disorder averaging. The de- 
tails of the approach for the dilution-only problem are given 
in Ref. 6 and results for M(x)/M(Q) are shown in Fig. Q] 
The modification of this method for the model concerns 
changes in the p- and d-wave scatterings off the impurities, 
while the s-wave contribution can be shown to be unaffected 
by the frustrating terms fisft . The advantage of this method 
is that both the x- and J' Zn , J Zn -dependence of the order pa- 
rameter can be studied systematically. 

QMC simulations were performed using the stochastic se- 
ries expansion method ll 811 to find the staggered magnetiza- 



(5) 



The frustrating interactions generally present a serious diffi- 
culty due to the negative sign problem [19], which becomes 
more serious (exponentially) with increasing number of frus- 
trating bonds and inverse temperature j3 = J/T. Here we 
focus on the system with a single impurity, where the sign 
problem is completely local (independent of the lattice size) 
and, in the range of frustration of interest here, manageable 
down to sufficiently low temperatures to draw conclusions 
about ground states of relatively large lattices. For all lattices 
considered below, Ml is well converged already at j3 = 16, 
and in the following we use this value. We assume the usual 
size-dependence of the staggered magnetization J^T 



21H: 
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As a check, we calculated M for a pure AF using the same 
system sizes and (3 = 16, which gave M w 0.3066 for the 
extrapolated value, only slightly below high-precision result 
obtained using much larger lattices; M — 0.30743(1) ll22ll . 

We are interested in the suppression of M(x) at low impu- 
rity concentrations; M(x)/M(0) w 1 - xR(0), where R(0) 
is the slope of the M vs x curve at x — > 0. We introduce a 
finite-size analog of the slope function as 



Rl = 



1 



1 



1 



Ml{1) 

M L (0) 1-x 



(7) 



where L is the lattice size, M^(l) is the staggered magnetiza- 
tion of the system doped by one impurity, and x = 1 /L 2 for 
that case. The normalization by (1 — x) is necessary to convert 
M found from (0 to weighting M relative to the amount of 
magnetic sites J8|]. Fig.|4fa) shows the size dependence of the 
slope in both the unfrustrated J-only model and the frustrated 
model with J Zn = 2J' Zn = 0.07J and 0.1J. The latter set is 
still below the frustration beyond which the sign problem be- 
comes too serious. As expected, the frustration increases the 
slope substantially. Given the non-linearity of the data, there 
is some uncertainty in the l/L—>0 extrapolation. We here use 
a linear fit to L > 8 data; see Fig.|4ja). Due to some remain- 
ing nonlinearities (which are seen clearly for sizes L < 8), the 
extrapolation may slightly under-estimate the slope i?(0), and 
should therefore be considered a lower bound. 

To obtain a direct quantitative measure of the effect of frus- 
tration we consider the difference ARl — R L — R L of the 
slopes from (O for the model (01 with (R* ) and without (R°) 
frustrating terms, respectively, for lattice sizes L and several 
J' Zn and J Zn . We perform finite-size extrapolations to the 
thermodynamic limit as above for each set of couplings; ex- 
amples are shown in the inset of Fig. |4|a). An unexpected 
finding is that two J Zn bonds suppress the order at almost 
the same rate as four J' Zn bonds of the same strength, as ev- 
idenced by both the QMC and the T-matrix results. Com- 
bining that with the systematically larger values of J'L from 
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FIG. 4: (Color online) (a) The slope Rl, Eq. 10 vs 1/L in systems 
with and without frustration (,9= 16), solid symbols are extrapolated 
values. Inset: slope increase ARl due to frustration, (b) Same as 
inset in Fig. [TJ with the T-matrix results for J'% n = 2J^ n = 0.07 J 
and QMC extrapolations of R(0) from Figfia) for jg n = 2J' Zn = 0, 
0.07 J, and 0.1J. QMC data provide lower bounds of the slopes. 



luted AFs and other doped Mott insulators. The character 
of the percolation transition should change as a result of the 
frustrating interactions across the impurities. The impurity- 
doping of spin chains should introduce weaker links between 
the broken pieces. Recent experiments on doped frustrated 
J1-J2 systems B25I1 should be affected by the same mecha- 
nism. The proposed impurity-induced frustrating interactions 
should persist at the finite doping and may produce a pair- 
breaking mechanism in the doped CUO2 planes. 
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the three-band model calculations shows that this interaction 
is particularly important. 

Fig-Htb) shows that the experimental slope R(0) « 1.1 is 
matched by the T-matrix results already at J Zn = 2J' Zn = 
0.07 J (J'tot = 0.28 J) where we fixed the ratio to 2 according 
to the three-band model results. QMC results from Fig. Ufa), 
extrapolated to x — ► 0, are shown in Fig. |4|b) for the same and 
larger values of frustration, J Zn = 2J' Zn =0.1 J ( J' tot =0.4 J). 
They yield the lower bounds for the slopes i?(0) ~0.7 and 0.9, 
respectively I23I1 . The results suggest that a stronger frustra- 
tion, close to the latter data set or somewhat larger, is present 
in the real Zn-doped CuC»2 plane. This is still a reasonably 
modest amount of frustration, well within the window sug- 
gested by the three-band model calculations. 

Using the same QMC analysis, we have also investigated an 
alternative mechanism of the enhanced order suppression by 
impurities, in which the strength of the J bonds in the vicinity 
of Zn is reduced by the lattice distortion [24]. We found that 
changing the bond strength by 15% (SJtot — 1.2J) changes 
the slope R(0) by at most a few percent, ruling out the lattice- 
distortion mechanism of the order suppression as a viable al- 
ternative to our theory. 

Conclusions. — We have proposed that impurity-doped 
strongly correlated systems develop significant frustrating in- 
teractions, which are absent or negligible in the corresponding 
undoped system, due to electronic degrees of freedom at the 
scale less than the Hubbard-C/. Applying this mechanism to 
the problem of the doped non-frustrated Heisenberg model, 
relevant to Zn-doped cuprates, we have found this effect to be 
the key to resolving discrepancies between experiments and 
earlier theories. Our analytical and numerical results agree 
well with each other. We estimate the total amount of frustra- 
tion to be > 0.4 J per impurity in the Zn-doped CUO2 plane. 
In light of our quantitative theoretical results, further high- 
precision experiments at low doping are called for. 

Outlook.-Om theory has far-reaching consequences for di- 
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